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1. Introduction 

Supersymmetry (SUSY) is a leading candidate 
for a theory of physics beyond the Standard 
Model. However, it is clear that if SUSY exists, it 
must be broken. The mechanism by which SUSY 
is broken is yet unknown and there is a growing 
list of possible scenarios. In these scenarios, the 
SUSY spectrum is described by a handful of pa- 
rameters, generally defined at the SUSY break- 
ing scale; RGE running of sparticle masses and 
coupling constants yields predictions for the mass 
spectra and decay patterns of the various sparti- 
cles at energy scales relevant for colliders or cos- 
mology. However, these SUSY breaking scenarios 
are restrictive and predict specific phenomenolo- 
gies for colliders and cosmology that may not rep- 
resent the full range of possible SUSY signatures. 

Here, we study the MSSM more broadly with- 
out assumptions at the high scale. We re- 
strict ourselves to the CP-conserving MSSM (i.e., 
no new phases) with minimal fiavor violation 
(MFV)[lJ. Additionally, we require that the first 
two generations of sfermions be degenerate as mo- 
tivated by constraints from fiavor physics. We 
are then left with 19 independent, real, weak- 
scale, SUSY Lagrangian parameters: the gaug- 
ino masses Mi ^2,3, the Higgsino mixing param- 
eter /i, the ratio of the Higgs vevs tan/3, the 
mass of the pseudoscalar Higgs boson m^, and 
the 10 squared masses of the sfermions (m^i^a. 



^ni,37^Ji 37^fi a^and mei,3). We include inde- 
pendent A-terms only for the third generation 
(^5, At^ and A^-) due to the small Yukawa cou- 
plings for the first two generations. This set of 19 
parameters has been called the phenomenological 
MSSM (pMSSM)[2j. 

To study the pMSSM, we performed a scan 
over this 19-dimensional parameter space as- 
suming fiat priors for the specified ranges [3j: 
100 GeV < < 1 TeV; 50GeV < |Mi,2,/i| < 
ITeV; 100 GeV < M3 < 1 TeV; \At^t,r\ < 
ITeV; 1 < tan/3 < 50; 43.5 GeV < tua < 1 TeV. 
We randomly generated 10^ points in this param- 
eter space and subjected them to an exhaustive 
set of existing theoretical and experimental con- 
straints. We also performed a similar scan with 
log priors (with slightly different mass ranges) to 
gauge the infiuence of priors on our results and 
found these to be negiglible[3j. We then gener- 
ated SUSY spectra utilizing SuSpect2.34[2]. 

2. Theoretical and Experimental Con- 
straints 

We now discuss the theoretical and experimen- 
tal constraints that we applied to the generated 
parameter space points; for more details, one 
should consult I3I4I . 



* Presented at the Dark Matter Conference, 9-11 Feb 2009, 
Arcetri, Florence, Italy. SLAC-PUB-13731. Work sup- 
ported in part by the Department of Energy, Contract 
DE-AC02-76SF00515. 



2.1. Theoretical Constraints 

We demanded that the sparticle spectrum not 
have tachyons or color or charge breaking min- 
ima in the scalar potential and also required that 
the Higgs potential be bounded from below with 
consistent electroweak symmetry breaking. We 
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assume that the LSP, which wih be absolutely 
stable, be a conventional thermal relic and iden- 
tify the LSP as the lightest neutralino. 

2.2. Low Energy Constraints 

The code micrOMEGAs2.20[5j was used to 
evaluate the following observables for each point 
in the parameter space: Ap, the decay rates for 
b ^ and Bs and the ^ — 2 of the 

muon. In addition, we evaluated the branching 
fraction for B ^ tv following [6j and [7J. We al- 
lowed a large range for the SUSY contribution 
to ^ — 2 due to the evolving discrepancy between 
theory and experiment [8j. We implemented con- 
straints from meson- antimeson mixing [9j by as- 
suming MFV[T], imposing first and second gener- 
ation mass degeneracy, and demanding that the 
ratio of first /second to third generation squark 
soft breaking masses differ from unity by no more 
than a factor of 5. 

2.3. Accelerator Constraints 

LEP data at the Z pole shows that charged 
sparticles with masses below Mz/'^ are unlikely. 
This also holds for the lightest neutral Higgs bo- 
son. Data from LEPIILIOJ indicates that there are 
no new stable charged particles of any kind with 
masses below 100 GeV. We also require that any 
new contributions to the invisible width of the Z 
boson be < 2 MeV[llJ. 

Results from sparticle searches at LEPII pos- 
sess numerous caveats. We implement a lower 
limit of 92 GeV on first and second genera- 
tion squark masses [12 and 95 GeV on the sbot- 
tom mass [13 , provided that the gluino is more 
massive than the squarks and the mass differ- 
ence (Am) between the squark and the LSP is 
> 10 GeV. We demand that the lightest stop 
mass be greater than 95(97) GeV[14] if the stop 
can (cannot) decay into Wbxi- The right-handed 
sleptons must have masses greater than 100, 95, 
or 90 GeV for selectrons, smuons, and staus re- 
spectively, as long as the condition 0.97msiepton > 
mLSP is satisfied. These bounds are also ap- 
plied to left-handed sleptons, when the neutralino 
t— channel diagram may be neglected in the case 
of selectrons. Chargino masses be greater than 
103(95) GeV, provided that the LSP-chargino 



mass splitting is Am > (<)2 GeV[14j. If 
the lightest chargino is dominant ly Wino, this 
limit only applies when the electron sneutrino 
t— channel diagram is negligible. The LEP Higgs 
Working Group[15j, provides five sets of con- 
straints on the MSSM Higgs sector, which are es- 
sentially limits on the Higgs- Z coupling times the 
Higgs branching fraction for various final states. 
We employ SUSY-HIT [16j to analyze these. In 
addition, we included a theoretical uncertainty on 
the calculated mass of the lightest Higgs boson 
of approximately 3 GeV [17] when applying these 
constraints. 

We also employ constraints from searches at 
the Tevatron. Restrictions on the squark and 
gluino sectors arise from the null DO multijet plus 
missing energy search [18]. We generalize their 
analysis, rendering it model independent, by gen- 
erating multijet plus missing energy events for 
our models using PYTHIA6.4[19j as interfaced to 
PGS4 [20j which provides a fast detector simula- 
tion. We weigh our results with K factors com- 
puted using PROSPINO2.0 [21]. Analogously, 
we employ constraints from the CDF search for 
trileptons plus missing energy [22], which we also 
generalize to the fuh pMSSM. D0[23j has ob- 
tained lower limits on the mass of heavy stable 
charged particles. We take this constraint to be 
m^+ > 206|/7i^p + m\Uih? GeV at 95% CL 
for charginos, where the matrix entries Uiw and 
Uih determine the Wino/Higgsino content of the 
lightest chargino. CDF and DO also have analyses 
that search for light stops and sbottoms[24j; these 
searches are difficult to implement in a model- 
independent pMSSM context and thus we exclude 
models with light (m < m^) stops or sbottoms. 

2.4. Astrophysical Constraints 

There are two constraints from considering 
the LSP as a long-lived relic. As noted above, 
we demand that the LSP be the lightest neu- 
tralino. We also require, following the 5 year 
WMAP measurement [25j of the relic density, that 
^^^|lsp ^ 0.121. In not employing a lower bound 
on 1^/i^|lsp for our models, we acknowledge the 
possibility that even within the MSSM and the 
thermal relic framework, dark matter may have 
multiple components. However, in discussing re- 
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suits below, we will also examine a subset of mod- 
els for which 0.100 < nh'^\LSP < 0.121. 

We also obtain constraints from direct dark 
matter searches [26 . Generally, the strongest con- 
straints come from the spin-independent WIMP- 
nucleon cross sections, hence we only imple- 
ment bounds on our models from these; inspec- 
tion of the spin-dependent WIMP-nucleon cross 
sections in our models confirms that this ap- 
proach is reasonable. Both spin-independent and 
spin-dependent cross sections were calculated us- 
ing micrOMEGAs2.21[5j. We implement cross 
section limits from XENON10[27], CDMS[28j, 
CRESST I[29J and DAMA[30J data. It should 
be noted that many of our models predict a 
value 1^/i^|lsp which is less than that observed 
by WMAP and supernova searches. We thus scale 
our cross sections to take this into account. 

3. Results 

As noted above, we randomly generated 10^ 
parameter space points {i.e., models) in a 19- 
dimensional pMSSM parameter space using flat 
priors. Only ~ 68.5 • 10^ of these models satisfy 
all the constraints listed above. The properties of 
these models are described in much greater detail 
in [3] . Here we will discuss the attributes of these 
models which are most important astrophysically. 

Figure [1] presents a histogram of the masses of 
the four neutralino and two chargino species in 
our models. The lightest neutralino is, of course, 
the LSP. The LSP mass lies between 100 and 250 
GeV in over 70% of our models. Generally mod- 
els with a mostly Higgsino or Wino LSP have a 
chargino with nearly the same mass as the LSP; as 
sufficiently light charginos would normally have 
been detected at LEP or the Tevatron, there are 
fewer models with such LSPs with mass ^ 100 
GeV. 

The identity of the nLSP is shown in Figure O 
The lightest chargino is the nLSP in about 78% 
of the models; this is due to many models having 
Wino or Higgsino LSPs, and the generally small 
mass splitting between a mostly Wino or Higgsino 
neutralino and the corresponding chargino. The 
second lightest neutralino is the nLSP ~ 6% of 
the time. These will generally be models with a 
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Figure 1. Distribution of neutralino masses (top 
panel) and chargino masses (bottom panel) for 
our set of models. 



dominantly Higgsino LSP. Note also that while 
neutralinos or charginos are the nLSP in the vast 
majority of cases, there are 10 other sparticles 
each of which is the nLSP in > 1% of our models. 
Scenarios in which these sparticles are the nLSP 
may lead to interesting signatures at the LHC [3T] . 

Figure [3] displays the LSP mass value as a func- 
tion of the LSP-nLSP mass splitting. Am, our 
models for each identity of the LSP. It is in- 
teresting that these models have a smaller Am 
than is often considered; 80% of our models have 
Am < 10 GeV, 27% have Am < 1 GeV, and 3% 
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Figure 2. Number of models in which the nLSP 
is the given sparticle. 



tence of such a region shows the difficulty of mak- 
ing model independent statements about sparticle 
masses or other SUSY observables. 

We have seen that within our model set the 
nLSP can be almost any SUSY particle and the 
corresponding Am can be small for these cases. 
Thus specific models in our set describe qualita- 
tively most of the conventional long-lived spar- 
ticle scenarios. Long-lived stops or staus (as in 
GMSB), gluinos (as in Split SUSY) as weh as 
charginos (as in AMSB) all occur in our sam- 
ple. We also have long-lived neutralinos, as does 
GMSB, however these are the X2 case. In 

addition to models which, to some extent, cor- 
respond to these well-studied scenarios, we also 
have models with long-lived selectrons, sneutri- 
nos and sbottoms. 



have Am < 10 MeV. As one can see from Fig- 
ure O this occurs largely, but not exclusively, in 
models with a chargino nLSP. This is again due 
to the many models where the LSP is nearly pure 
Wino or Higgsino. 

There are a number of interesting features in 
this figure. The mostly empty square region 
which appears on the lower left-hand side of Fig- 
ure [3] is due to the fact that models with chargino 
nLSPs in this mass and Am range have been ex- 
cluded by the Tevatron stable chargino search. 
Non-chargino nLSPs are not eliminated by this 
search {e.g., the production cross section for slep- 
tons in this range is too small to be excluded by 
the Tevatron search). It is perhaps worth noting 
that a stable heavy charged particle search at the 
LHC, corresponding to those done at the Teva- 
tron, would be able to exclude or discover the 
models with heavier chargino nLSPs and small 
values of Am (corresponding to ~ 12% of our 
model set). 

Another interesting feature in this figure is 
the bulge for 0.1 GeV < Am < 2 GeV and 
^LSP ^ 100 GeV. This region exists because 
these values of Am are large enough that at LEP 
or the Tevatron, the produced chargino would 
decay in the detector, but the resulting charged 
tracks would be too soft to be observed. The exis- 
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Figure 3. Mass splitting between nLSP and LSP 
versus LSP mass. The identity of the nLSP is 
shown as well. (The LSP is always the lightest 
neutralino in our set of models). 



Figures [H and [5] display the gauge eigenstate 
content of the LSPs in our model set. We note 
that most LSPs are relatively pure eigenstates, 
with models where the LSP is Higgsino or mostly 
Higgsino being by far the most common. About 
one quarter of our models have Wino or mostly 
Wino LSPs, while just over one-sixth have Bino 
or mostly Bino LSPs. Within mSUGRA, the LSP 
is, in general, nearly purely Bino; this suggests 
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Figure 4. The distribution of LSP gaugino eigen- 
state types as a function of the LSP mass (top 
panel) or the LSP-nLSP mass difference (bottom 
panel). Note that each LSP corresponds to three 
points on this figure, one each for its Bino, Wino, 
and Higgsino fraction. 



that most of our models are substantially differ- 
ent from mSUGRA. We note that one would ex- 
pect the LSP be a pure eigenstate fairly often in 
a random scan of Lagrangian parameters, since if 
the differences between Mi,M2, and /i are large 
compared to M^, then the eigenstates of the mix- 
ing matrix will be essentially pure gaugino and 
Higgsino states. 




bino fraction 

Figure 5. Wino/Higgsino/Bino content of the 
LSP in the case of fiat priors. Note that, as else- 
where in the paper, |^i2p, and |^i3p + 
|Zi4p, where Zij is the neutralino mixing ma- 
trix, give the Bino, Wino, and Higgsino fractions 
respectively. 
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3.1. Relic Density 

We did not demand that the LSP, in any given 
model, account for ah of the dark matter, rather 
we required only that the LSP relic density not 
be too large to be consistent with WMAP. More 
specifically, we employed 1^/i^|lsp < 0.121. Fig- 
ure [6] shows the distribution of 1^/i^|lsp values 
predicted by our model set. Note that this dis- 
tribution is peaked at small values of ^h^\hSP- 
In particular, the mean value for this quantity in 
our models is ^ 0.012. We note that the range 
of possible values of 1^/i^|lsp is found to be much 
larger than those obtained by analyses of specific 
SUSY breaking scenarios [32j. We display the pre- 
dictions for 1^/i^|lsp versus the LSP mass and ver- 
sus the nLSP - LSP mass splitting in Figure [71 
Figure [3 makes it clear that ^h^lhSP generally 
increases with the LSP mass, but a large range 
of values for the relic density are possible at any 
given LSP mass. The empty region in Figure [3 
where 1^/i^|lsp ~ 0.001-0.1 and mLSP ^ 50-100 
is due to the fact that, in general, LSPs which 
are mostly Higgsino or Wino give lower values 
of ^/i-^ I LSP 5 and, as noted above, there are fewer 
Higgsino or Wino LSPs in this mass range. This 
figure also shows that small mass differences can 
lead to large dark matter annihilation rates. 
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Figure 7. Distribution of 1^/i^|lsp as a function of 
the LSP mass (top panel) or the LSP-nLSP mass 
splitting (bottom panel. 
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Figure 6. Distribution of ^h^\hSP for our models. 



3.2. Direct Detection of Dark Matter 

As noted above, we calculate the spin- 
dependent and spin-independent WIMP-nucleon 
cross sections using micrOMEGAs 2.21 [5 . These 
data give the possible signatures in our model set 
for experiments that search for WIMPs directly. 
As these experiments measure the product of 
WIMP-nucleon cross sections with the local relic 
density, the cross section data presented in the fig- 
ures below are scaled by ^ = fth'^ \ lsp /^h'^\ wm ap • 
To date, these experiments generally provide a 
more significant bound on the spin-independent 
cross section, and hence we will focus on those. 

Figure [8] presents the distribution for the 
scaled WIMP-proton spin-independent cross sec- 
tion versus relic density for our model sample. 
As one would expect, larger values of the cross 
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section are generally found at larger values of 
^^^Ilsp- However, even for relic densities close 
to the WMAP value, £,(Jp^si is seen to vary by al- 
most eight orders of magnitude. These ranges for 
i^p,si are much larger than those from mSUGRA 
as calculated, e.g.^ in [33] . 



WIMP-Proton Ggi Versus Relic Density 
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Figure 8. Distribution of scaled WIMP-proton 
spin-independent cross section versus the LSP 
contribution to relic density for our models. 



Figure [9] shows the scaled WIMP-proton spin- 
dependent and spin-independent cross sections as 
a function of the LSP mass. The constraints 
from XENON10[27J and CDMS[28] are also dis- 
played. As noted above, to take the uncertain- 
ties in the theoretical calculations of the WIMP- 
nucleon cross section into account, we allowed for 
a factor of 4 uncertainty in the calculation of the 
WIMP-nucleon cross section. Table 3 in Ref.[4j 
gives the fraction of models that would be ex- 
cluded if the combined CDMS/XENONIO cross 
section limit were improved by an overall scaling 
factor. Note that our inclusion of the theoreti- 
cal uncertainties does not significantly modify the 
size of our model sample. 

We find that the range of values obtained for 
these cross sections covers the entire region in 
cross section/ LSP space that is anticipated from 
different types of Beyond the Standard Model 
theories in the above reference. This suggests 




100 200 300 400 500 SOO 700 



LSP Mass Versus WIMP-Proton Ogi 
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Figure 9. Distributions of scaled WIMP- 
proton spin-dependent cross section and spin- 
independent cross sections versus LSP mass in 
our models. In the spin-independent panel, the 
CDMS and XenonlO bounds are shown. 
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that we cannot use direct detection experiments 
to distinguish between e.g. SUSY versus Little 
Higgs versus Universal Extra Dimensions dark 
matter candidates in the absence of other data. 

In Figure [TOl we compare the WIMP-proton 
and WIMP-neutron cross sections in the spin- 
dependent and spin-independent cases. The spin- 
independent cross sections are seen to be fairly 
isospin independent; this is not the case, however, 
for the spin-dependent cross sections. 
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Figure 10. Comparison of the WIMP-neutron 
and WIMP-proton cross sections. The spin- 
dependent (independent) cross sections are shown 
in the top(bottom) panel. 



3.3. Indirect Detection of Dark Matter 

The PAMELA collaboration has recently 
claimed an excess in the ratio of cosmic ray 
positrons to electrons observed at energies ^ 10 
GeV[34j. Here we employ DarkSUSY 5.0.4[35] 
to calculate this ratio for our model sample and 
compare these results with the PAMELA data. 

In general, for a thermal relic dark matter can- 
didate to reproduce the PAMELA data, its signal 
rate must be multiplied by a boost factor [36j. In 
nature, such a boost factor could result from, e.g.^ 
a local overdensity. The boost factor in that case 
would be the square of the ratio between the den- 
sity of dark matter in the region from which one 
is sensitive to cosmic ray positrons and electrons 
to the universe as a whole. 

In our analysis, we use four propagation mod- 
els which are present in darkSUSY: the model of 
Baltz and Edsjo[37], that of Kamionkowski and 
Turner [38j, that of Moskalenko and Strong [39j, 
as weh as GALPROP|40j. These are referred 
to in the following figures as "BE","KT","MS", 
and "GAL", respectively. Interestingly, we find 
that the extent to which the positron/electron 
fiux ratio predicted by our models matches the 
PAMELA data can be highly sensitive to the 
choice of propagation model parameters. We will 
explore this further in future work[41J. 

The differential positron flux as a function of 
energy for a random sample of 500 models from 
our set are shown in Figure [Til Here we assume a 
boost factor of 1; the normalization of the curves 
takes into account the fact that for many of these 
models 1^/i^|lsp < ^wmap- 

We next determine how well the predicted 
positron fluxes for these models agree with the 
PAMELA data, allowing for the possibility of a 
boost factor. To do this, we flnd the value for 
the boost factor (with the restriction that it be 
< 2000) which minimizes the fo^ the flt of each 
model's prediction to the PAMELA data. In cal- 
culating the we consider only the seven high- 
est energy bins, as at lower energies solar modu- 
lation is expected to play a major role[34j. Fig- 
ure [12] shows the and corresponding boost fac- 
tor for these 500 random models. Note that there 
are four data points for each model, so there are 
actually 2000 data points in this flgure. We then 
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Figure 11. Expected flux spectrum of positrons 
from neutralino annihilation in the halo for 500 
randomly selected models. For each model, there 
are three curves, one for each of three propagation 
models (as shown in the legend and defined in the 
text). The dotted black line is the expected back- 
ground of positrons from non-SUSY processes. 
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Figure 12. The distribution of P^r degree 
of freedom versus the choice of boost factor that 
minimized this quantity for 500 randomly se- 
lected pMSSM models in our model set. These 
quantities have been determined for each of four 
propagation models. Only boost factors less than 
2000 were considered; this explains the large num- 
ber of models for which the x^-minimizing boost 
was 2000. 




Figure 13. Positron/ electron flux ratio versus 
energy for the pMSSM models for which the 
per degree of freedom with the x^-i^iaximizing 
boost was less than 10.0 for three of the four prop- 
agation models. Curves are shown for all four 
propagation models. 



display the positron to electron flux ratio, for the 
models with a low value of x^, as a function of en- 
ergy in Figure [131 and note the reasonable agree- 
ment with the data for some models. 

Since the flux from WIMP annihilation scales 
as (^^/i-^lLSp/^^^lwiviAp)^, we might expect to 
improve the match to the PAMELA data using 
models from our sample for which I^/i^Ilsp ^ 
^^^IwMAP- To test this, we examine the pre- 
dicted positron flux for 500 random models with 
^^^Ilsp > 0.100; these fluxes are shown in Fig- 
ure [Ml with no boost factor. We then again flnd 
the boost factor that minimizes the x^ of the 
positron to electron flux ratios with respect to 
the seven highest energy PAMELA bins; these are 
shown in Figure [151 Here, we note that there are 
many more models for which the x^-iiiinimizing 
value for the boost factor is < 2000 and there 
are many more points for which the x^ value is 
low. The positron to electron flux ratios for these 
models, including the boost factor, are shown in 
Figure [161 

It appears that some of our models do a rea- 
sonably good job of fltting the PAMELA positron 
data, especially in the case where 1^/i^|lsp lies 
fairly close to the WMAP value. For most mod- 
els, describing the PAMELA data requires large 
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Figure 14. Expected flux spectrum of positrons 
from neutralino annihilation in the halo for 500 
randomly selected models for which I^/i^Iwmap > 
^^^|lsp > 0.10. For each model, there are three 
curves, one for each of three propagation models 
(as shown in the legend). The dotted black line is 
the expected background of positrons from non- 
SUSY processes. 
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Figure 15. The distribution of per de- 
gree of freedom versus the choice of boost factor 
that minimized this quantity for 500 randomly 
selected pMSSM models with 1^/i^|wmap > 
^^^|lsp > 0.10. These quantities have been de- 
termined for each of four propagation models. 
Only boost factors less than 2000 were consid- 
ered; this explains the large number of models 
for which the x^-minimizing boost was 2000. 




Figure 16. Positron/ electron flux ratio versus 
energy curves for pMSSM models for which the 
per degree of freedom with the x^-maximizing 
boost was less than 5.0 for three of the four prop- 
agation models. Curves are shown for all four 
propagation models. 



boost factors, however this is also a fairly generic 
feature of attempts to explain PAMELA and 
ATIC data in terms of WIMP annihilation [36]. 
There are however, many models which give rel- 
atively low x^ per degree of freedom in the fit 
to the data with relatively small boost factors. 
We will study this further in future work [41j. A 
study of the corresponding predictions for the the 
cosmic ray anti-proton flux is also underway. 

4. Conclusions 

We have generated a large set of points in 
parameter space (which we call "models") for 
the 19-parameter CP-conserving pMSSM, where 
MFV has been assumed. We subjected these 
models to numerous experimental and theoreti- 
cal constraints to obtain a set of ~ 68 K models 
which are consistent with existing data. We at- 
tempted to be somewhat conservative in our im- 
plementation of these constraints; in particular 
we only demanded that the relic density of the 
LSP not be greater than the measured value of 
Vtb? for non-baryonic dark matter, rather than 
assuming that the LSP must account for the en- 
tire observed relic density. 

Examining the properties of the neutralinos in 
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these models, we find that many are relatively 
pure gauge eigenstates with Higgsinos being the 
most common, followed by Winos. The rela- 
tive prevalence of Higgsino and Wino LSPs leads 
many of our models to have a chargino as nLSP, 
often with a relatively small mass splitting be- 
tween this nLSP and the LSP; this has impor- 
tant consequences in both collider and astropar- 
ticle phenomenology. 

We find that, in general, the LSP in our models 
provides a relatively small (~ 4%) contribution 
to the dark matter, however there is a long tail 
to this distribution and a substantial number of 
models for which the LSP makes up all or most 
of the dark matter. Typically these neutralinos 
are mostly Binos. 

Examining the signatures of our models in di- 
rect and indirect dark matter detection experi- 
ments, we find a wide range of signatures for both 
cases. In particular, we find a much larger range 
of WIMP-nucleon cross sections than is found in 
any particular model of SUSY-breaking. As these 
cross sections also enter the regions of parameter 
space suggested by non-SUSY models, it appears 
that the discovery of WIMPs in direct detection 
experiments might not be sufficient to determine 
the correct model of the underlying physics. As a 
first look at the signatures of these models in indi- 
rect detection experiments, we examined whether 
our models could explain the PAMELA excess in 
the positron to electron ratio at high energies. We 
find that there are models which fit the PAMELA 
data rather well, and some of these have signif- 
icantly smaller boost factors than generally as- 
sumed for a thermal relic. 

The study of the pMSSM presents exciting new 
possibilities for SUSY phenomenology. The next 
few years will hopefully see important discover- 
ies both in colliders and in satellite or ground- 
based astrophysical experiments. It is important 
that we follow the data and not our existing prej- 
udices; hopefully this sort of relatively model- 
independent approach to collider and astrophysi- 
cal phenomenology can be useful in this regard. 
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